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Abstract. We study the resonant decay of the primordial Standard Model Higgs conden¬ 
sate after inflation into SU(2) gauge bosons on the lattice. We find that the non-Abelian 
interactions between the gauge bosons quickly extend the momentum distribution towards 
high values, efficiently destroying the condensate after the onset of backreaction. For the 
inflationary scale H = 10® GeV, we find that 90% of the Higgs condensate has decayed after 
n ~ 10 oscillation cycles. This differs significantly from the Abelian case where, given the 
same coupling strengths, most of the condensate would persist after the resonance. 
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1 Introduction 

After the detection of Standard Model (SM) Higgs [1, 2] there has been a growing interest 
in the cosmological role of the Higgs field 4>. Intriguingly, the Higgs might even act as the 
inflaton if its potential at high energies is dominated by a coupling to the spacetime curvature 
of the form [3]. The connection to the low energy SM regime of the scenario is, 

however, not unique due to its non-renormalizability [4-6] and the required large coupling 
^ > 10^ could be subject to questions of naturalness [7-9]. 

However, the Higgs could also be of cosmological interest in several other ways. If the 
Higgs potential remains close to the SM form up to the inflationary scale R*, during inflation 
the Higgs would be a light spectator field, as discussed in [10-17]. In such a case inflation 
needs to be driven by some new physics beyond SM. Like for all light fields during inflation, 
both the Higgs mean field as well as the local field value will be subject to fluctuations. 
The latter are isocurvature perturbations during inflation while fluctuations of the mean 
field generate an effective Higgs condensate. The condensate survives inflation and sets 
specific non-vacuum initial conditions for the hot Big Bang epoch [10, 12]. The typical 
magnitude of the Higgs mean field is, provided inflation lasts long enough, of the order of 
. Depending on details of physics beyond SM the primordial Higgs condensate could have 
significant observational impacts ranging from imprints in GMB to baryogenesis and dark 
matter [11, 18-21]. 

Theoretical self-consistency of the setup requires a stable vacuum. As is well known, 
the pure SM vacuum in flat space becomes unstable above energies fic r-u 10^*^ GeV for the 
measured best fit values of the strong coupling constant, the top mass, and the Higgs mass 
M/j ~ 125 GeV [10, 14-17, 22-26]. However, as discussed in [10, 27], during inflation the 
stability depends crucially also on the Higgs-curvature coupling This coupling is 

necessarily generated by radiative corrections even if ^ would be set to zero at some scale. A 
one-loop investigation [27] shows that vacuum stability can be maintained even within the 
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SM up to the maximal inflationary scale -fT* < 10^^ GeV consistent with the tensor bound 
[28, 29], provided that ^ew ^0.1 at the electroweak scale. 

Thus quite generically, immediately after inflation the Universe features a primordial 
Higgs condensate. In the case of Higgs inflation, in which the Higgs itself sets the dynamics, 
the Universe is dominated by the Higgs condensate, and its decay time directly determines 
the reheating temperature [30, 31]. On the other hand, within the SM, and in a broad class 
of its extensions, the condensate contributes very little to the energy density of the Universe. 
At the onset of hot Big Bang its large displacement from the SM vacuum could nevertheless 
have important ramifications for the subsequent evolution of the Universe [11, 18, 20, 21]. 
Thus even if the Higgs were a mere spectator during inflation, it can affect post-inflationary 
physics until the condensate eventually decays. In both cases a detailed understanding of the 
Higgs condensate decay is crucial in order to properly address the observational imprints. 

At zero temperature the dominant decay channel for both the SM Higgs and the non- 
minimally coupled Higgs inflaton is the non-perturbative production of gauge bosons through 
a parametric resonance [12, 30-33]. The investigation of the resonance is complicated by the 
non-Abelian dynamics of the gauge field which have not been carefully explored so far. In the 
regime of broad resonance, realized for SM running of the gauge couplings, the non-Abelian 
terms are small during the first stages of the resonance but rapidly grow important as the 
exponential production of gauge fields starts to backreact on the dynamics of the Higgs 
condensate. The efficiency and duration of the resonant Higgs decay therefore cannot be 
reliably estimated without properly accounting for the non-Abelian couplings. In the limit 
of a narrow resonance, which could occur in specific extensions of the SM, the non-Abelian 
features affect the resonance dynamics already before the dynamical backreaction of gauge 
fields [32]. 

In this work we investigate the resonant decay of the Higgs into gauge bosons using 
numerical lattice simulations and accounting for the full non-Abelian dynamics. In particular, 
we compute the decay time of the subdominant SM Higgs condensate generated during 
inflation. 

The paper is organized as follows. In Section 2 we present the framework to investigate 
the resonant Higgs decay and make some analytical estimates. In Section 3 we present the 
results of our full lattice computation of the Higgs decay which constitutes the main part of 
the paper. Finally we conclude in Section 4. 


2 Resonant Higgs decay in Standard Model 


During and shortly after inflation the energies are significantly higher than the electroweak 
scale. Thus the Standard Model Higgs potential can be approximated by the quartic term 
alone. The relevant part of the SM action for our discussion is 


5 = - y d^x + A($tck)(<kt(i))2] | (2.1) 


with the kinetic terms 


F;, = , 

G^iU = ^ flBy — , 

= fv. - igAy^ - CD . 


( 2 . 2 ) 
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Here $ is the Higgs doublet and and are the SU(2) and U(l) gauge fields, respectively. 

We assume that the inflationary stage is driven by some physics beyond the SM. To 
reheat the universe, the inflaton field, or fields, should arguably be coupled to the SM degrees 
of freedom through non-gravitational interactions. We assume these couplings are small and 
use the pure SM running for the Higgs and gauge couplings. More importantly, we will 
be investigating the decay of the Higgs field assuming there is no thermal background. If 
the inflaton has already decayed at this stage, we are therefore implicitly assuming that 
the interaction rate between its decay products and the SM fields is small compared to the 
Hubble time over our epoch of interest. On the other hand, one could argue that the decay 
rate of the inflaton field would typically be much smaller than the decay rate of the Higgs 
condensate. In such a case the oscillating inflaton field would give rise to an effectively matter 
dominated universe. 


2.1 Initial conditions from inflation 

Given the SM potential, the SM Higgs is both a light field and energetically subdominant 
during inflation [11, 12]. Its super-horizon dynamics can be investigated using the stochastic 
formalism [34, 35] . Starting from a generic field configuration at some point during inflation, 
the Higgs distribution relaxes to the equilibrium given by P{h) ~ exp(—87r^l//3ff^) within 
a time scale of N^ei. — 100 e-folds [36]. Provided that inflation lasted somewhat longer we 
can therefore quite generally assume that the value of the Higgs condensate in the observable 
universe is to be drawn from the equilibrium distribution. 

After the end of inflation the typical mean Higgs field value over the observable universe 
is thus given by 

h* = . (2.3) 

The inflationary stage thus generates a primordial Higgs condensate /i*, setting specific non¬ 
equilibrium initial conditions for the Hot Big Bang epoch. Here we are implicitly assuming a 
low inflationary scale PI* < 10® GeV such that the Higgs coupling to curvature scalar 
can be neglected [27]. 


2.2 Onset of the resonant Higgs decay 

The Higgs condensate becomes effectively massive as = 3A*/i* and it starts to oscillate 
around the vacuum expectation value h = 0. In terms of the rescaled field x = a/i and 
conformal time dr = a~^dt the equation of motion reads 


X + Ax^--X = 0. (2.4) 

a 

The last term can in general be neglected. In a radiation dominated background it is iden¬ 
tically zero and for matter domination it becomes negligible soon after the onset of Higgs 
oscillations. The expansion of space can therefore be scaled out from the dynamics and the 
problem essentially reduces to flat space. 

The resonant production resulting from coherent oscillations of scalar fields has been 
extensively studied in the literature [37-40]. Neglecting self-interactions, the weak transverse 
components of the gauge fields obey the Lame differential equation 


d^A'^ 


2 , 2 
K -|- (7 cn 



= 0 . 


(2.5) 
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1 fr for the W bosons 
AxLc [ ^ 4 / for the Z boson 

and q is the resonance parameter characterizing the strength of particle production. For the 
pure Standard Model this parameter is always greater than unity, which puts the system 
in the broad resonance regime [32]. The Lame equation exhibits instability regions where 
the solutions are amplified exponentially resulting in explosive production of gauge boson 
quanta. 

The exponential production of gauge fields could in principle be washed out by rapid 
perturbative decays of the generated gauge fields. This is indeed what happens in the early 
stages of reheating after the Higgs inflation [30, 31]. For the energetically subdominant SM 
Higgs this is however not the case. The perturbative decay rate of W bosons into any pair 
of fermions is small compared to the characteristic time scale of Higgs oscillations 

— = \/3327r^ = 0(0.01) . (2.7) 

V A 

Similar estimates hold for the Z bosons. The perturbative decays are therefore irrelevant 
compared to the resonant production of gauge fields and can be neglected in what follows. 

2.3 Backreaction and limitations of the analytical approach 

In the discussion of the previous section we have ignored the non-Abelian interactions of the 
gauge fields as well as the backreaction of the produced particles on the dynamics of the Higgs 
field. This is justified in the beginning of the resonance where the gauge field occupation 
numbers are small. 

As the occupation numbers grow larger than one, the non-Abelian terms backreaction 
effects eventually start to dominate the dynamics. Using the Hartree approximation one can 
estimate the onset of backreaction by evaluating the time when [32] 

2 2/2 

~ 3Ax2 . (2.8) 

The backreaction destroys the coherence of the Higgs oscillations and eventually shuts down 
the resonance [38-40] . The dynamics in this regime is entirely dominated by interactions and 
its analysis requires a full non-linear lattice computation. Before turning to the numerical 
results in the next section let us however present some qualitative comments on the role of 
non-Abelian interactions. 

For Abelian fields the distribution of occupation numbers would retain its form peaked 
around k ~ g^^^/2 and essentially a vacuum configuration for k ^ The low energetic 
gauge bosons and Higgs particles k < would be inefficient in fragmenting the remnants of 
the Higgs condensate. Consequently, it would seem that a significant part of the condensate 
would survive long after the backreaction. 

The non-Abelian interactions of the W bosons however completely change the picture. 
As the occupation numbers around the resonant peak k ~ {‘1 grow much larger than 

unity the W scatterings and annihilations into Z bosons rapidly generate excitations with 
higher momenta. This generates a bath of gauge bosons extending to which would 

appear to efficiently destroy the remnants of the Higgs condensate, in sharp contrast with 






the Abelian case. As we will discuss below, this behaviour is exactly what we see in our 
numerical computation. 

We should note here that also annihilations of the gauge bosons into two fermions 
participate in the Higgs decay. Indeed, the vacuum cross sections of the processes WW ff 
and WW —>■ WW both scale as a ~ and the production of light fermions is also 

kinematically favoured, as discussed in the context of preheating after the Higgs inflation 
[30, 31]. However, at the hnal stages of the resonance the non-Abelian interactions WW —?■ 
WW^WW —)■ Z are significantly favoured by the large occupation numbers nw,z >1- In 
our lattice computation we will therefore neglect the decay channels into light fermions. 


3 Full lattice computation of the Higgs decay 


Motivated by the considerations above, we simulate the SM Higgs dynamics on lattice during 
the period over which exponential production of gauge bosons is taking place. As the Higgs 
decay proceed dominantly through the production of non-Abelian W bosons we simplify 
the problem by investigating a SU(2) gauge-Higgs system instead of the full SU(2) x U(l) 
symmetry. Technical details of the lattice simulation can be found in Appendix A. The 
techniques we use are similar to those employed in previous numerical studies of preheating 
involving gauge fields [41-43]. 

Note that the system under consideration is conformal and so the effect of expansion 
can be neglected. Therefore there is no mass scale, other than that set by the initial value 
of the Higgs field. 

Our initial conditions for the gauge field are Gaussian white noise with strength 

Q)A]{k\ 0)) = - k') . (3.1) 


rather than the quantum vacuum. To avoid needing to carry out projection to satisfy the 
Gauss law initially, we assume that the gauge field conjugate momentum is initially zero. The 
gauge field equilibrates quickly, in any case, and loses all memory of these initial conditions 
as the particle number starts to grow. We then take as z = 0 the time when the particle 
occupation number at the peak of the resonance exceeds = 0.1. 

We have confirmed that the normalization ^ of the gauge field initial conditions only 
has a logarithmic impact on the time taken for the system to equilibrate. For the purposes 
of the current paper, then, the above Gaussian noise initial conditions are adequate. 

The lattice simulations are carried out in temporal gauge, which is equivalent to setting 
Aq = 0, however we must further fix lattice Coulomb gauge to recover particle numbers. 
When we wish to measure particle numbers in the gauge field we use standard gauge fixing 
techniques (see Appendix A.l) and then measure the connected two point functions [43-45] 


(>l“(k.().4‘{-k,i)) 

{£f(k,t)£‘(-k,i)) 


^ab 
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'^)D^(k,t) + '^DE(k,t) 


With Df^{k,t) = 0 in Coulomb gauge we define the gauge field particle number as 


(3.2) 

(3.3) 


n^{t) = ^JD^{k,t)D§{k,t). 


( 3 . 4 ) 











This quantity is not gauge invariant and so we also present the power spectrum of the energy 
density in the gauge field in our results. This observable does not show the clear sharp 
resonant peak but it does allow one to identify backreaction and the onset of equilibration. 

For comparison we also plot the results of a ‘non-compact’ U(l) simulation with equiv¬ 
alent initial conditions^. A comprehensive numerical study of the resonant production of 
Abelian gauge bosons was recently carried out in Ref. [33], and so we will not go into detail 
here. 

In general we use lattices of volume 128^. With gauge fixing, each simulation at this 
volume typically required four hundred CPU-hours; in total the results of this paper required 
less than ten thousand CPU-hours. On the infrared side we varied the lattice volume to 
ensure that the parametric resonance was not unacceptably cut off at long wavelengths, and 
confirmed that the results of this paper do not depend on this. It is less straightforward to 
establish that the physics does not depend on rescaling the system, because it is conformal 
and there is no mass scale. 

To provide a numerical estimate of the time taken for the Higgs field to backreact, we 
define the backreaction time z^r as the time when the amplitude of the Higgs oscillations 
falls to 90% of its initial value fi*. We emphasise, however, that this time is logarithmically 
dependent on the amplitude ^ of the initial conditions for the gauge field. The simulations 
presented here are however more than sufficient to demonstrate the correct qualitative be¬ 
haviour. 

After the backreaction on the Higgs commences, the system starts to equilibrate and 
the equation of state for the scalar field approaches the equipartition value (1/3). We also 
see indications of thermalisation. However, because we do not have a mass scale in our 
simulations, we cannot give a definitive estimation of the thermalisation time. This must be 
left to future work that fully incorporates expansion. 

3.1 Results 

To give a broad resonance peak within our lattice volume we rescaled the system so that the 
lattice spacing is unity and work with rescaled lattice Higgs field $ = 2%. 

In Figure 1 we show how the backreaction time varies with q for the noise amplitude. 
For the Standard Model changes of q = 5^/(4A) correspond to changes of the inflationary 
scale which determines the values of the running coupling constants, identifying the 
renormalization scale as ^ = /i* ~ H*; see Eq. (2.3). For example, taking = 10® GeV 
yields with the SM running and best fit inputs the coupling values g ~ 0.6 and A = 0.15 at 
the end of inflation [23]. This gives the resonance parameter q = g'^/(4A) ~ 6. The peaks 
(slowest to backreact) are associated with the values of q where the first (lowest-momentum) 
resonant band encountered is farthest from k = 0 . 

Figure 2 depicts the time-evolution of the non-Abelian gauge field particle numbers. 
For the sake of comparison, we have also shown the time evolution in an Abelian case with 
equivalent Higgs and gauge coupling strengths. The initial evolution is similar in both cases 
but a drastic difference is observed as the particle numbers become much larger than unity. 
After this point the non-Abelian distribution is rapidly extended towards higher momenta 
through mutual interactions of the gauge bosons. In the Abelian case, where the interactions 

^Note that a lattice simulation of ‘compact’ U(l) would not correctly show the onset of backreaction as the 
photon field Ai{x) would be symmetric under transformations Ai{x) —>■ Ai(x) + 2nn. For the same reason, 
it is not sufficient to only initialise one colour of the SU(2) gauge field to get an ‘abelian’ comparison to the 
nonabelian case. 
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Figure 1: Plot of the backreaction time Zbr against q for the non-Abelian theory. It seems that 
some fine tuning is required to substantially change the Higgs decay time. Also shown is the smallest 
K that lies within a resonance band at a given q (for a full contour plot of resonance bands see e.g. 
Ref. [3:^1' 



Figure 2: Comparison of particle numbers for abelian (dashed lines) and non-Abelian (solid 
lines) simulations (g « 6.12) evaluated at evenly-spaced time intervals. Equilibration in the abelian 
case is much slower and particle numbers remain many orders of magnitude smaller. 


are absent the distribution on the other hand remains strongly peaked around the momenta 
K < excited by the resonance. Similar behaviour can also be observed in the gauge field 
energy densities illustrated both for the Abelian and non-Abelian case in Figure 3 

The bath of non-Abelian gauge particles with momenta extending to k rapidly 

destroys the Higgs condensate through scattering processes. This can be observed in Figure 4 
which depicts the evolution of the various energy components of the Higgs and SU(2) gauge 
helds. Comparing Figures 2 and 4 we clearly see that the Higgs condensate rapidly starts 
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(a) 2 = 8.4 (b) z = 16.4 (c) z = 24.4 (d) z = 32.4 


Figure 3: Slices of the gauge field energy density, comparison between Abelian and non-Abelian 
simulations, q = 6.12, at successive time-slices. Backreaction for the non-Abelian case has clearly 
started before z « 24 and the non-Abelian self-interaction means that by z « 32 the field configuration 
is no longer ‘smooth’. In contrast, the Abelian simulation is still dominated by resonant structures at 
z « 32. The factor of three difference in the legend normalizations is due to the additional non-Abelian 
colour degrees of freedom. 



Figure 4: Time evolution of the individual energy components for a typical non-Abelian simulation 
(g«6.12). 


to decay after the gauge boson particle numbers at the resonance peak have ceased their 
exponential growth. The resonance shuts off and the remnants of the Higgs condensate 
decay completely within a few oscillation cycles. 
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Figure 5: Comparison of scalar potential energies for Abelian and non-Abelian simulations (g « 
6.12). The decay of the Higgs condensate is substantially faster in the non-Abelian case. 




Figure 6: Plots showing (at left) the particle numbers n„ and (at right) the power spectrum of 
the non-Abelian gauge field, at equal times (intervals of 5z = 4) throughout a typical simulation 
(g « 6.12). The power spectrum closest to z = 0, at z = 0.4 is shown in both plots as a red dashed 
line. 


In the Abelian case, a significant fraction 0(1} of the condensate remains even after 
the end of resonance. This is illustrated in Figure 5 where we compare the evolution of 
the condensate energy in the Abelian and non-Abelian case. As discussed above, in the 
Abelian case the gauge field distribution remains peaked around the long-wavelength modes 
k < excited by the resonance. Correspondingly, the gauge boson scatterings of the Higgs 
field are not sufficient to destroy the condensate and the Higgs remains close to its initial 
non-equilibrium configuration. 

Finally, in Figure 6 we show the development of particle number for a non-Abelian 
simulation at finer time intervals than is possible in figure 2, as well as the gauge field power 
spectrum for the same intervals, to give a gauge-invariant picture of the energy deposited in 
the bosons. 
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Note that the results presented here (and in particular Figure 2) bear a strong resem¬ 
blance to Figure 7 in Ref. [46]. The underlying physics of equilibration and, presumably, 
subsequent thermalisation is very similar. Corresponding results have also been obtained for 
tachyonic preheating in Ref. [43]. 

4 Conclusions 

In this work we have studied the non-perturbative decay of the primordial Higgs condensate 
after inflation into SU(2) gauge bosons which are produced through parametric resonance 
once the condensate leaves slow-roll and starts to oscillate. We performed a full lattice 
simulation of the resonant production of gauge bosons focusing in particular on the detailed 
dynamics after the onset of backreaction and carefully accounting for the effects arising from 
non-Abelian interactions. 

Resonant decay of the condensate results in the production of particles within particular 
bands of momenta in the infrared region. However, we have found that the scattering and 
decay processes resulting from non-Abelian interactions rapidly extend the momentum dis¬ 
tribution into the ultraviolet, which efficiently destroys the remnants of the Higgs condensate 
after the onset of backreaction. This is in sharp contrast to the purely Abelian case where 
the distribution remains in the infrared for much longer and the produced particles do not 
efficiently destroy the condensate so that a sizeable part of it survives after the resonance is 
terminated. 

The time for backreaction to set in and the Higgs condensate to decay is found to be 
logarithmically dependent on the inflationary scale 77* which determines the Higgs initial 
conditions. For 77* = 10® GeV we find that 90% of the Higgs condensate has decayed after 
n ~ 10 oscillation cycles. 

After the decay is completed the system starts to evolve towards equilibrium. Energy 
is rapidly transferred from the initially occupied IR modes towards higher momenta and the 
system appears to evolve towards a stationary state. This resembles the self-similar behaviour 
observed in earlier lattice studies of preheating [47, 48] and also studies of heavy ion collisions, 
see e.g. [49, 50]. The evolution of this stage towards the full thermal equilibrium can however 
not be analysed within a classical simulation which fails in describing the UV modes. We 
leave a more detailed investigation of the thermalisation for a future work. 

We have emphasized that unless the SM Higgs sector is significantly modified by new 
physics, a Higgs condensate is inevitably generated during inflation. Thus it is important to 
understand the details of its decay as it may be significant for subsequent physics. The Higgs 
condensate sets specific out-of-equilibrium initial conditions for the hot Big Bang epoch which 
could have significant ramifications ranging from primordial perturbations and baryogenesis 
to non-thermal production of dark matter [11, 18-21]. In all these cases the decay time 
of the Higgs condensate is a crucial factor affecting the observational impacts. Our results 
indicate that, in the absence of thermal bath, once the Higgs condensate unfreezes and starts 
to oscillate it decays completely after 0(10) oscillation cycles. 

While we have considered the decay of the condensate within the context of the Standard 
Model it would be interesting to extend the analysis to the case of Higgs inflation, where the 
condensate plays the role of the inflaton through a large non-minimal coupling to gravity. 
While order of magnitude corrections to previous estimates [30, 31] are not expected, a proper 
inclusion of non-Abelian interactions appears crucial to precisely determine the predicted 
reheating temperature of the scenario. Other extensions of the Standard Model may also 
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be of interest. For instance, in some extensions the resonance might turn out to be narrow 
rather than broad. In such a case the non-Abelian nature of the interactions is expected to 
become important already before the onset of backreaction [32] . 
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A Details of lattice implementation 

The lattice Hamiltonian is 



+ '^^\24)\x)(l){x) — 2Re(l)\x)Ui{x)4>{x + i) ^ U(</>^(x)(/>(x)) (A.l) 


X 


where vr and P are the conjugate momenta for the Higgs and gauge helds, Ui{x) are the link 
variables corresponding to the gauge field (see below) and (j) is the Higgs held. The equations 
of motion are then 


(/>(t + 6t, x) = (j){t, x) + 6t7r{t + St/2, x) 

Ui{t + St, x) = exp —i^P/‘{x,t + St/2)cr‘^St Ui{t,x) 
7r{t + St/2, x) = 7r(t — St/2, x) 


(A.2) 

(A.3) 


-I- St ^Ui{t, x)(f){t, x + i) — 2(f>{t, x) + Uj{t, X — i)4>{t, x — i) 


_ i 


dV 

5(/>t 

(A.4) 


P/^{t + St/2,x) = P/^{t — St/2,x) + St ^Re 4 /'{t,x + k)Uj,{t,x)i(T'^<j){t,x) 



+ia'^Uk{t,x)ul{t,x + k - i)ul{t,x - i)Ui{t,x - i) 


(A.5) 


on the lattice, with Gauss law 


G{x) = Tr i a^Pi{t + St/2, x) — i cr^uj{t, x — i)Pi{t + St/2, x — i)Ui{t, x — i) 


i 


-|-2g(Re TT^{t + St/2, x)ia^(j){t, x) =0. (A.6) 


- 11 - 



















See for example Refs. [41-43] for more details. 

Our measurements of non-Abelian gauge boson particle number use a lattice approxi¬ 
mation to Af{x,t). Given the parallel transporter is defined through 


Ui{x, t) = exp 


a 




we deduce 

A°‘(x,t) = -Tv a°‘Ui{x,t) 
9 

to second order in the lattice spacing. 


(A.7) 


(A.8) 


A.l Gauge fixing in SU(2) 

Our equations of motion are in temporal gauge, and so while at some point the system may 
be in Coulomb gauge, the time evolution of the system will generally take it away from 
this form. We hx Coulomb gauge on the lattice [51, 52], when we wish to measure particle 
numbers. To do this we must find a lattice gauge transformation g{x) G SU(2) such that 

F<^[U] = :^Y.^eTvU!ix) (A.9) 

x,i 

is stationary, where 

Uf{x) = g{x)Ui{x)g\x + i). (A.IO) 

We look for a maximum of this expression by carrying out a combination of simulated an¬ 
nealing and overrelaxation steps, a strategy that is widely used for lattice quantum field 
theory simulations where fixing to Coulomb, Landau or Abelian gauge is required (see for 
example Refs. [52-54]). We perform 10000 sweeps of simulated annealing combined with over¬ 
relaxation. The divergence A(x) of the equivalent of the continuum gauge field, Eq. (A.8), 
computed through 

- ?)) (A.ll) 

i 

can be normed to give a measure 9 of the quality of the gauge fixing, 


After our gauge fixing procedure we typically obtain 0 ~ 10 given an initial value around 
10 “ 2 . 

The above algorithm as implemented in our code is fully parallelised so the need to fix 
the gauge to produce meaningful particle number results is, overall, not a significant burden. 

To obtain the electric held Ef{x,t) note that the Coulomb gauge condition is time- 
independent. We therefore carry out the gauge hxing procedure - with the same g{x) - on 
plaquettes at two sequential timesteps U^{x, t) and Uf{x, t+5t) and use these to approximate 
the electric held Ei{x,t) through 


exp 




Uf{x,t + 6t)Uf'^{x,t). 


(A.13) 
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Note that the gauge transformation g{x), while time independent, does not necessarily bring 
configurations on any other timestep into Coulomb gauge, so our estimation of the gauge-fixed 
electric field could have 0 {6t) errors. 

Having obtained lattice approximations to and Af{x,t) we can compute Eqs. 

(3.2-3.3) and obtain particle numbers for the system at a given time from Eq. (3.4). 

Note that, as demonstrated previously in Ref. [43], at large particle numbers the gauge 
fixing procedure does not substantially affect the final results. One could infer the same 
conclusions as we reach in this paper from the transverse ‘particle number’ obtained without 
gauge fixing; the discrepancy is around 10%. 

A.2 Gauge fixing in U(l) 

The same principle is applied here, except that a simple gradient flow procedure easily deter¬ 
mines the transformation Ai{x,t) + A(x,t) that satisfies the condition diAi = 0. 

The same gauge transformation A can then be applied to the following timestep to obtain 
Ei{x,t) in a manner analogous to that given above. 

A.2.1 Numerical tests 

We confirmed that our backreaction times and other principal results were logarithmic in the 
initial energy density We also checked that the results did not depend on the timestep 
size or the lattice volume V (we tested 128^ and 256^, both of which capture the parametric 
resonance well). While we can test the large-volume and infrared robustness of our results, 
the system is conformal and there is no mass scale. Therefore we cannot take the continuum 
limit or draw any conclusions about physics once particle numbers in the UV start to grow 
after backreaction. 

The simulations should always satisfy the Gauss law. The worst-case total Gauss law 
violation G{x)‘^/V in the non-Abelian case was around 10“^°, and occurred around the 
time of backreaction. 
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